function [X,Y] = distribution(mir_wide)
%采用满排列的方式 构造不同镜子形状的镜子分布

dR=mir_wide+5;
R=100:dR:350;%计算每一圈的半径
N=length(R);%圈数
diff_angle=(mir_wide+5)./R;

X=[];
Y=[];

for i=1:N
    if mod(i,2)~=0%奇数圈
        [Xn,Yn]=pol2cart(pi/2:-diff_angle(i):-pi/2,R(i));%半圈的镜子分布
        X=[X,Xn];
        Y=[Y,Yn];
    elseif mod(i,2)==0%偶数圈
        [Xn,Yn]=pol2cart(pi/2-diff_angle(i-1)/2:-diff_angle(i):-pi/2,R(i));%半圈的镜子分布
        X=[X,Xn];
        Y=[Y,Yn];
    end
end
X=X';
Y=Y';

end